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ABSTRACT 

Kepler-78b is one of a growing sample of planets similar, in composition and size, to the 
Earth. It was first detected with NASA’s Kepler spacecraft and then characterised in more 
detail using radial velocity follow-up observations. Not only is its size very similar to that 
of the Earth (1.2f?©), it also has a very similar density (5.6 g cm -2 ). What makes this planet 
particularly interesting is that it orbits its host star every 8.5 hours, giving it an orbital distance 
of only 0.0089 au. What we investigate here is whether or not such a planet could have been 
perturbed into this orbit by an outer companion on an inclined orbit. In this scenario, the outer 
perturber causes the inner orbit to undergo Kozai-Lidov cycles which, if the periapse comes 
sufficiently close to the host star, can then lead to the planet being tidally circularised into a 
close orbit. We find that this process can indeed produce such very-close-in planets within the 
age of the host star (~ 600 — 900 Myr), but it is more likely to find such ultra-short-period 
planets around slightly older stars (> 1 Gyr). However, given the size of the Kepler sample and 
the likely binarity, our results suggest that Kepler-78b may indeed have been perturbed into 
its current orbit by an outer stellar companion The likelihood of this happening, however, is 
low enough that other processes - such as planet-planet scattering - could also be responsible. 
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1 INTRODUCTION 

Analysis of data from NA SA’s Kepler spacecraft l lBorucki et al.1 
12010 : iBatalha et all 12013l) indicates t hat planets with radii simi¬ 
lar to that of the Earth are common JPetigura, Marcv & Howard 
l2013l : lDressing & Charbonneaij|2013l ). Recently it was announced 
that one of the Kepler targets (Kepler-78) showed a 0.02% decline 
in brightness that was associated with a plane t with a radius of only 
1.16 ± 0.19.R® iSanchis- Oieda et all |201 3h- Follow-up observa¬ 
tions, using HARPS-N l iCosentino et al . 2012) and the High Reso¬ 
lution Echelle Spectrometer (HIRES) Ivogt et alll994l ). confirmed 
that this is indeed a planet with a mass of about 1.86M® and a d en- 
sity of about 5.6 g cm -3 jPene et alj 20131 ; [Howard et al] 20131) . 

Of course it is fascinating that we are now detecting plan¬ 
ets with sizes and densities similar to that of the Earth, but what 
makes this planet particularly interesting is that it has an orbital pe¬ 
riod of only 8.5 hours, meaning that it is orbiting at a distance of 
only 0.0089 au from its parent star. Quite how such a planet can 
end up in such an orbit is very uncertain. It almost certainly could 
not have formed where it now resides, as the temperature in the 
disc in that region would hav e been too high even for dust grains 
to condense aBell et alJI 1997D . It could potentially have migrated 
inwards through disc migration. However, such low-mass planets 
would migrate in the gapless, Type I regime l l Ward 179971) which 


* E-mail: wkmr@roe.ac.uk 


is typically thought to be so fast l lTanaka. Taleuchi & WardfeoQTl : 
IKlev & Cridal 12008 f) that it would seem unlikely that such objects 
could be left st randed so close to their pa rent stars. P opulati on syn¬ 
thesis models llda & Linll2008l : Mordasini, Alibert & Benj l20091) 
typically assume a reduced Type I migration rate. 

Alternativel y, suc h close- in planets could be scattered onto 
eccentric orbits dFord & Rasid 120061) that ar e then circularise d 
through tidal interactions with the parent star dRasio et~aDl 19961 ). 
It has indeed been suggested that if such a process were to oper- 
ate, we sh ould be seeing som e very short period hot super-Earths 
dSchlaufman. Lin & IdauOlOh . so this could be an explanation for 
the origin of Kepler-78b. However, even this study suggested that 
typical orbital periods would be greater than the 8.5 hour orbital 
period of Kepler-78b. 

Another mechanism for forming close-in planets, related to 
dynamical interactions in multi-planet systems, is for the planet to 
undergo Kozai-Lido v cycles driven by a stell ar companion on a 
highly inclined orbit dKozail[7962l : iLidovll 1962l) . If the eccentricity 
is sufficiently large, so that the periastron becomes very small, the 
planet’s o rbit may be circularised through tidal interactions with its 
host star dWu & Murravll2003l : lFabrvckv & Tremaind[2007l) . What 
we want to investigate in this paper is whether or not this process 
could indeed explain the origin of Ke pler-78b. Given that binarity 
amo n gst solar -lik e star s is quite high iDuquennov & Mavonl 19911 : 
lAbt & Willmarthl2006h it seems likely that this could play a role in 
producing close-in, Earth-sized planets. 
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2 Ken Rice 


In this paper we present results from a series of Monte Carlo 
simulations in which we consider how a planet with a mass and 
density the same as that of Kepler-78b, but initially orbiting be¬ 
tween 0.5 and 2 au, is influenced by perturbations from a binary 
stellar companion. We also include the influence of tides, which 
would allow the orbit to circularise if the eccentricity becomes suf¬ 
ficiently large, and the influence of general relativistic and apsidal 
precession. The paper is organised as follows : in Section 2 we 
present equations of motion, in Section 3 we describe the basic 
setup of the problem, in Section 4 we discuss the results, and in 
section 5 we discuss the results and draw some conclusions. 


stellar wind. FromlBarker & Qgilviel <120091') we have 
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2 EQUATIONS OF MOTION 


The goal is to evolve an inner binary (planet and star) under the 
influence of tidal interactions, perturbing accelerations from stel¬ 
lar and planetary distortions due to tides and rotation, perturbing 
accelerations from a third body, and general relativistic apsidal 
precession. To quadrupole order, these equations were first pre- 
sented by (jE ggleton & Kiselev a||200l|) and can also b e found in 
IWu & Murray! < ?2003h and Fabrvckv & Tremain~el d 20071) . 

_ Rather than using the equations in Eggletonj&Kixelevi] 

d200lh. we ’ ve im plemented those from iBarker & Qgilvie J2009h 
and I Barker] d201 ll) which are regular at e = 0. We want to evolve 
an inner system (planet + host star) where the bodies have masses 
M s and M p , radii R s and R p , and in which the orbit has an 
eccentricity e, semi-major axis a, and orbital angular frequency 
n — \JG(M a + M p )/a 3 . The vector quantities that we want to 
evolve are, therefore, the spin of the parent star fi s , the spin of the 
planet f2 p , the eccentricity of the inner orbit e, and angular mo¬ 
mentum vector of the inner orbit h = r x r = na 2 y/ 1 — e 2 h. 

We’ve built our model by considering, initially, only the equa¬ 
tions that evolve these quantities through tidal dissipation and a 


where I s and I p are the moments of inertia of the star and planet, 
S\wmd represents the stellar wind, and ^ = M S M P /(M S + M p ) 
is the reduced mass of the inner system. We also need to define the 
tidal friction timescales for the star and planet (tf s and t f p ), which 
depend on the star and planet’s tidal quality factors ( Q' s and Q ' p ), 
and the functions of the eccentricity. 
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We also include the contributions due to an additional outer body 
(b) of mass M 0 , orbital angular frequency n 0 , semi-major axis 
a 0 , and eccentricity e 0 . Additionally, we add contributions from 
quadrupolar distortions of the inner star and planet due to tidal and 
rotational bulges (qs and qp), and we include general relativistic 
apsi dal precessio n ( GR ). The equations, shown below, are taken 
from lBarkej d201 il l 
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Cb M s + M p + M 0 n 4(1 - e 2 ) 1 / 2 (l - e 2 ) 3 / 2 ' (19) 

In Equation 118b . k s and k p are the inner star and planet’s tidal love 
numbers. 


2.1 Octupole terms 

It now appears that expanding the equations only to q uadrupole or¬ 
der jnay not be ap propria te for many systems jNaoz et all 1201 it 
[Naoz. Farr & Rasiq 12012h . so we’ve also included the octupole 
terms. This allows us to consider situations in which the outer 
body’s mass is comparable to that of the inner planet, and to con¬ 
sider situations in which the outer orbit is eccentric. 

We’re unable to write the octupole terms in a way that is reg¬ 
ular a t e = 0, so have implemented the form in jMardling & Linl 
l2002lt . The octupole contributions are 
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where n is a unit vector that is perpendicular to the plane of the 
outer body’s orbit (not to be confused with n and n D , which are the 
angular frequencies of the inner and outer orbits), and 
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where R is the co-ordinate of the outer body, and the co-ordinate 
frame is defined by the basis vectors ( e,q,h ), with q = h x e. The 
other unit vectors above are 7?i - R ■ e, 7+ = R- q, R 3 = R- h. 
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4 Ken Rice 


2.2 Integrating the outer orbit 

The octupole terms described above need the co-ordinate of the 
outer body, which we determine by solving for the eccentric 
anomaly, E. This can be done by iterating the following equations 
until dE is below a threshold (we use 10 -12 ) 

-{E - e 0 smE - l) 

ah, = -—- 

1 — e Q cos E (22) 

E = E + dE, 

where l is the mean anomaly [l = n 0 (t — P)] , with P the orbital 
period and t the time since the completion of the last full orbit of 
the outer body. In all of our simulations, we fix the outer body to 
lie in the xy plane and so its co-ordinates are then 

Rx = a 0 {cosE — e 0 ) 

R y = do\f\--e I sin E (23) 

Rz = 0, 

where a 0 and e 0 are the semi-major axis and eccentricity of the 
outer orbit. In this work, we neglect perturbations on the outer orbit. 


2.3 Stellar wind 

Without a stellar wind, or with a very weak stellar wind, 
it is possible that tidal interactions between the star and 
planet can result in the plane t being trapped in a close orbit 
dDobbs-Dixon. Lin & Mardlinal 120041 ). Most stars, however, have 
winds that continue to remove angular momentum, and so once 
tidal interactions become significant, we would typically expect the 
planet to continue spiralling in towards the central star. We imple- 
ment here a very simple magnetic braking form for the stellar wind 
i Weber & David 1 19671 : iKawalet! 1 19881 : ICollier Cameron & Jiankel 
119941) so that in the unsaturated regime, the stellar wind term is 

^swind = V 3 S , (24) 

where k w is the braking efficiency coefficient. In the saturated 
regime (when Q. s > T2) this becomes 

f l.S'wi Tlfl n;n. (25) 

The braking efficiency coefficient, k w , is set so that the stellar rota¬ 
tion period matches that expected for the star being considered, and 
Cl is set to be 14flg. The vector associated with the stellar wind is 
always set so as to point in to opposite direction to that of the spin 
of the planet host star. 


2.4 Putting it together 

Ultimately we want to evolve the angular momentum, h, and ec¬ 
centricity, e, of the inner orbit, and the spins of the planet and its 
host star, T2 P and T2 S . The evolution of the stellar spin is determined 
by combining the stellar wind equations [Equations 1241 and 1251 
with Equation 0. The evolution of the spins of the star and planet 
[Equation 0] , both depend on the tidal evolution of the orbital 
angular momentum [Equation 0J. 

To evolve the angular momentum of the inner orbit, we need 
to add the contributions from tides [Equation (0}] , perturbations 
from an outer body expanded to quadupole and octupole order 
[Equations 1131 . and 1201 ] , and perturbations from distortions of 
the inner star and planet [Equation G3]. Similarly to evolve the 
eccentricity of the inner orbit, we combine the contributions from 
tides [Equation 0]. perturbations from an outer body expended 
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Figure 1. A figure showing the evolution of the semimajor axis (a - solid 
line) and the peri aps (a(l — e) — dash ed line) using initial conditions the 
same as thos e in|W u & Murray (2003). This was a code test that was also 
carried out bv lFabrvckv & Tre maine ( 20 071), and ou r results appear to match 
theirs. It doesn’t quite match I Wu & Murrayl 120031) but we can match their 
results if we ignore the term representing the apsidal precession due to the 
spin and tidal bulges of the planet. 


to quadrupole and octupole order [Equations l |l5K and 1211 , per- 
turbations from distortions of the inner star and planet [Equation 
1161 ] and general relativistic apsidal precession [Equation 1171 ]. 


2.5 Some basic tests 


Since this is a new code, we ran a few comparison tests to check 
that it was w orking properly. The first was that introduced by 
Wu & Murray! <20031) . It comprises a 1.1 Mq star with a 7.8Mj up 
planetary companion, the star having a radius of IRq and the 
planet having a radius the same as that of Jupiter. The initial stel¬ 
lar and planetary spin periods are, respectively, 20 days and 10 
hours. The inner system’s orbit has a semimajor axis of a = 5.0 au, 
and eccentricity e = 0.1, the tidal love numbers are k s = 0.028 
and k p = 0.51, and the tidal dissipation quality factors are Q' g = 
5.35 x 10' and Q' p = 5.88 x 10 5 . The system also has a l.lMg 
companion with a 0 = 1000 au, e 0 = 0.5 and with an orbital plane 
inclined at 85.6° to that of the plane of the inner orbit. 

Figure |T| shows the time evolution of the semi-major axis 
(dashed line) and periaps (solid line) of the system described a bove, 
and appears the same as that in lFabrvckv & Tremaine! <20071). who 
alsojterformed this test. It’s not quite the same as in Wu & Murray! 
<2003l) . but we can match their result if we remove the apsidal pre¬ 
cession due to the spin and tidal bulges of the planet (which can 
require very short timesteps and, hence, long integration times). 

The second test was primarily to check that the octupole terms 
had been properly implemented. In this test, taken from lNaoz et al.1 
<2013l) . we ignore the tidal evolution terms, the terms associated 
with the distortion of the inner star and planet due to their tidal 
bulges, and the effect of general relativistic apsidal precession. The 
system consists of an inner star of mass 1 Mq, a companion planet 
with mass lMj up , and an outer planet with mass lMg. The inner 
orbit has a semimajor axis of a = 6 au and eccentricity of e = 
0.001. The outer orbit has a semimajor axis of a 0 = 100 au, an 
eccentricity of e 0 = 0.6 and is inclined at 65° to the plane of the 
inner orbit. The argument of pericentre of the inner orbit is also set, 
initially, to 45° with the outer one set to zero. 

Figure [2] shows the time evolution of 1 — e for the system 
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Figure 2. The ev olution of 1 — e for a system with the same parameters 
as those used by iNaoz et ail J20 1 3l> and described here in the text. The 
solid line show s the evo lution when the octupole terms are included and it 
matches that of lNaoz et al.l J20 1 3t> . The dashed line shows how the system 
would evolve in the absence of the octupole terms. 


d escri bed above. The result appears identical to that in INaoz et alj 
d2013t ). The dotted line also shows how the system would evolve in 
the absence of the octupole terms. 


3 BASIC SETUP 

The system we want t o consider specifically is Kepler-78 
dSanchis-Oieda et aPl2013l l. The companion planet, with a mass of 
M p = 1.86AT® and radius of R p = 1.1737?®, is Earth-sized and 
has an Earth-like d ensity (p = 5.57 g cm~ 3 ) dPene et al.ll2013l : 
iHoward et al.ll2013h . The host star has a mass of M s = 0.81 M® 
and radius R s = 0.737M®, and the planet has an orbit that is cir¬ 
cular (e = 0) and orbits at a distance a = 0.0089 au (giving an 
orbital period of 8.5 hours). 

The system is thought to h ave an age between 600 and 900 
Myr dSanchis-Oieda et alj|2013r) . so here we run our initial simu¬ 
lations for 800 Myr. We set the stellar wind braking parameter to 
k w = 10 46 , which gives a stellar rotation speed of b etwee n 11 
and 1 2 days at t = 800 Myr, similar to that observed dPepe et al.1 
l2013h . A star with a mass similar to that of Kepler-78, however, 
does not spin down much in the first Gyr and so the stellar wind 
is probably not particular important here. The tidal love numbers 
are set to k s = 0.028 and k p = 0.51, and we consider stellar tidal 
quality factors of Q' s = 5 x 10 5 and Q' s = 5 x 10 6 . We’ll specify 
the tidal quality factor for the planet Q p , the planet’s initial orbital 
properties and the properties of the outer body, when we discuss 
the results of the simulations. 


4 RESULTS 

Since we want to consider if Kozai-Lidov cycles could explain 
the properties of Kepler-78b, our models are set up in the follow¬ 
ing way. We assume we have a planet with a mass and radius the 
same as that of Kepler-78b with an initial semi-major axis between 
a = 0.5 au and a = 2 au, with the semi-major axis chosen ran¬ 
domly in log a. The initial eccentricity is set to be e = 0.05, cho¬ 
sen because we’re assuming, here, that the planet has formed in a 


circumstellar disc in an almost circular orbit. We should acknowl¬ 
edge, however, that the initial eccentricity can have a signi ficant 
impact on the evolution of co-planar systems dLi et al J|20 1 4|) and, 
therefore, stress that our results only apply to systems in which the 
inner system has a low initial eccentricity. 

We also assume that there is an outer companion with a 
mass randomly chosen to be uniform between M 0 = 0.1 M® and 
M 0 = 1A7®, a semi-major axis chosen randomly in log a , between 
a 0 = 40 au and a 0 = 20000 au, and a randomly chosen eccentric¬ 
ity between e 0 = 0 and e 0 = 1. We then fix the outer companion's 
orbit to be in the xy plane and randomly ori entate the inn er orbit so 
that the mutual inclination, i, is isotropic l iWu. Murray & Ramsahail 
i2007b . We also randomly orientate the longitude of the planet’s as¬ 
cending node. By choosing such a hig h-mass c ompani on, we ’re es¬ 
sentially in the test particle regime dLithwick & Naozll201 ill . Such 
companions will also produce a large maximum ecce ntricity (for 
the inner orbit) than lower mass companions dTevssandier et al] 
l2013ll . As such, we mi ght expect a reasonably l a rge number of 
tidal disru ption events (INaoz. Farr & Rasiol l2012t iLi et al.l I20l4 
IPetrovichl uO 15). As such, our results only apply to a situation 

where the companion is of stellar mass. _ 

We also impose stability critera 1 Lithwick & Naozl 1201 ll: 


INaoz et al.l2013l:lMardling & Aarsethll200 ill and inisist that 


e G 


a 0 1- el 


< 0 . 1 , 


(26) 


and that 


— > 2.8 

a 



M 0 

M s +M p 


2/5 (1 + e 0 ) 2/s / 0.37 \ 

(1 - e D ) 6 / 5 V 180° ) ' 

(27) 


Equation ( 126b ensures that we are in the regime where the 
quadrupole and octupole terms dominate, while Equation ( 127b . in 
which i is the mutual inclination of the two orbits, ensures tha t 
the triple system is long-term stable ( iMardling & Aarsethil200ll ). 
Equation ( 127b is almost always satisifed for the initial conditions 
used here. 


4.1 Initial results 

The tidal quality factor for a terre strial planet is thought to lie be¬ 
tween Q p = 10 and Q' p = 500 dGoldreich & Soterll 1 9661) . Since 
we’re considering a young system in which the planet likely re¬ 
tains a lot of its initial internal heat, we assume a value at the 
top of this range (Q „ = 500), and als o a m ore ex treme case 
where Q' p = 5000 ( Henning. O’Connell & Sasselo\il2009l) . For 
the star, we assume tidal quality factors of Q' s = 5 x 10 s and 
Q' s = 5 x 10 6 , within the range ex pected for ex oplanet host stars 
fearaffe. Chabrier & Barman! I 2 OI 0 I : iBrown et al.l |201 lh . For each 
simulation we select the initial conditions as described above and 
evolve the system until t = 800 Myr, using a fourth-order Runge- 
Kutta integrator. We repeat this 10000 times for each set of param¬ 
eters, and the basic result is shown in Figure[3] The top panel is for 
Qp = 500 and the bottom for Q p = 5000. The solid line in each 
figure is for Q' s = 5 x 10 6 , the dashed line is for Q' s = 5 x 10 5 , 
and the vertical dash-dot line indicates the current semimajor axis 
of Kepler-78b. In each case, the number of planets still located be¬ 
tween 0.5 and 2 is very large and their distribution extends well 
above the limits shown on the y-axis. 

From Figure [3] it seems clear that it is possible for a planet 
to be perturbed into an orbit inside a = 0.01 au within 800 Myr. 
However, the numbers are typically small. For Q' p = 500 it is 10 
(Q' s = 5 x 10 6 ) and 4 (Q' s = 5 x 10 5 ), while for Q' p = 5000 it is 
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a (AU) 



a (AU) 


Figure 3. Histograms showing the final semi-major axes of the simulations 
with initial setup as described in the text and in which each simulations is 
stopped at t = 800 Myr. The top panel is for Q' p = 500 and the bottom 
is Q' p = 5000. The stellar tidal qualifty factors that we consider are Q' s = 
5 X 10® (solid line) and Q' s = 5 X 10 5 (dashed line). We consider 10000 
system in each case, with the planet starting with a between 0.5 au and 
2 au. After 800 Myr, there are between 10 and 4 planets surviving inside 
a = 0.01 au, depending on the values of Q' p and Q' s . The vertical dash-dot 
line indicates the current semimajor axis of Kepler-78b. In each case, the 
number of planets still located between 0.5 and 2 au is very large, and their 
distribution extends well above the limits shown on the y-axis. 


7 and 4 respectively. Even though the number of planets surviving 
inside a = 0.01 au is small, in most cases, a much larger number 
reach their Roche limit [a =_R E /0.462(M*/M P ) 1 ^ 3 = 0.0056 
au] J Faber. Rasio & Willemsll2005i ) and are assumed to be tidally 
disrupted and destroyed. With the exception of the Q' p = 500, 
Q’ s = 5 x 10® simulation (in which the numbers were small), in 
excess of 100 - out of a sample of 10000 - reached the Roche limit. 

Given that a large numbers of planets do become tidally de¬ 
stroyed, it is useful to know for how long a planet might exist in¬ 
side a = 0.01 au. Figure [4] shows four single planet simulations, 
one for each combination of Q' p and Q' s , each of which is run un¬ 
til the planet reaches its Roche limit. The amount of time such a 
planet spends inside a = 0.01 au depends, primarily, on the star’s 
tidal quality factor. For Q' s = 5 x 10®, the planet reaches the Roche 
limit in 480 Myr, while for Q' s = 5 x 10 5 it takes 48 Myr. There¬ 
fore, it would seem that for reasonable estimates of the star’s tidal 
quality factor, a planet such as Kepler-78b will only be detectable 
inside a = 0.01 au for a few hundred Myrs at most. 

Our initial results would therefore seem to suggest that it is 



Figure 4. A set of single planet simulations, each with a different Q' p or Q' a 
value. In each case we’ve chosen inital conditions that would indeed perturb 
the planet into a close orbit, and have run each simulation until the planet 
reaches its Roche limit (a = 0.0056 au in this case). The time a planet 
spends inside 0.01 au before reaching its Roche limit depends primarily on 
the star’s tidal quality factor and varies from 48 Myr ( Q' s = 5 X 10 s ) to 
480 Myr ( Q' s = 5 X 10®). 


possible for an outer companion to perturb a planet like Kepler-78b 
into a very close orbit (a < 0.01 au) within the age of the system 
(~ 800 Myr). However, the numbers are small, with at most 10 out 
of 10000 surviving inside a = 0.01 au at t = 800 Myr. 

4.2 Age of the system 

The previous simulations only considered the likelihood of a sys¬ 
tem with an age similar to that of Kepler-78, having a planetary 
companion with orbital properties similar to the of Kepler-78b. The 
results suggest it is possible, but probably rare. Additionally, the 
age d istribution of a sample of 950 Kepler object of interest host 
stars l lWalkowicz & Basril2013l) suggests that about 10% have ages 
less than 1 Gyr. This is also probably biased towards younger stars 
because of the way in which the sample was selected . Given that 
about 40% of these would have stellar companions d Rag ha van et al.1 
l20ld) would further reduce the likelihood of actually observing a 
Kepler-78b type system. 

To investigate how our results might depend on the age of the 
system, we reran our simulations with the same set of parameters as 
described above, but allowed the age of the star to vary, uniformly, 
from 500 Myr, to 2 Gyr. The resulting histograms are shown in 
Figure[5]and are very similar to those in Figure[3] There is a slight 
increase in the number surviving inside a = 0.01 au for Q' s = 
5 x 10®. In these runs there were 17 and 16 for Q' p = 500 and 
Q'p = 5000 respectively, compared to 10 and 7 when the age of the 
system was fixed at t = 800 Myr. For Q' a = 5 x 10 5 , the numbers 
are similar to the runs with the age fixed at t = 800 Myr. 

To further see the influence of the age of the system, we plot 
in Figure[6]the final semimajor axis of the planet against age of the 
system, for all those systems in which planets end up inside 0.05 
au. We only show, however, the results for Q' p = 5000, Q' s = 
5 x 10® as that produced the largest number of surviving planets 
inside a = 0.01 au. The first thing to note is that it appears more 
likely to detect such a planet for systems older than 1 Gyr (Kepler- 
78 has an age of between 600 and 900 Myr). However, Figure [6] 
does show 4 systems with an age < 1 Gyr, and with a planet inside 
a = 0.01 au. 
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Figure 5. Histograms showing the final semi-major axes for the simulations 
with initial setup as described in the text and in which the age of each system 
is randomly chosen to be between 500 Myr and 2 Gyr. The top panel is for 
Q' p = 500 and the bottom is Q' p = 5000. The stellar tidal qualifty factors 
that we consider are Q' s = 5 X 10 6 (solid line) and Q' s = 5 X 10 5 (dashed 
line). We consider 10000 system in each case, and the planets start with 
a between 0.5 and 2 au. In each case, the number of planets still located 
between 0.5 and 2 au is very large, and their distribution extends well above 
the limits shown on the y-axis. Depending on the values of Q p and Q' g , 
there are between 17 and 6 planets surviving inside a = 0.01 au. 


IWalkowicz & Basil i2013l) also suggest that maybe 20% of 
the Kepler tar gets ha ve ages less th an 2 Gyr. Kepler observed about 
150000 stars l lBorucki et all2010h . which suggests maybe as many 


as many as 30000 could have ages less than 2 Gyr. Candidates 
as small as Kepler-78b, however, are typically found around qui¬ 
eter - and therefore older - stars. Kepler is therefore incomplete 
for stars with high Combined Differential P hotometric Precision 
jBatalha et al.ll2013l : IChristiansen et~aT] 120 1 3l) and so the number 
of such planets is likely an underestimate. In the scenario shown 
in Figure [6] 17 - out of 10000 - survive inside a = 0.01 au. 
If we assum e t hat 40 % of those stars have stellar companions 
dRaghavan et al. 2010|), and that all of those stars could host ter- 
restrial planets Icassan et all2012l : lGreaves & Ricel2oi~i]) . then we 
might expect as many as 20 of the 30000 Kepler targets, with ages 
below 2 Gyr, to host such a ultra-close-in planet. Of course, our 
other simulations suggest that the number surviving could be as 
low as 2 (depending on the tidal properties of the star and planet), 
but given that the chance of such a system transiting is actually 
quite high (46%), observing such a system is still quite likely. 

Similarly, if we consider only those systems with ages below 


Figure 6. Figure showing the final semimajor axes for those planets that end 
up with a < 0.05 au, plotted against the age of the system. This simulations 
each considered 10000 systems in which the outer perturber was assumed 
to have a semimajor axis distribution that extended to a 0 = 20000 au. It’s 
clear that it is possible for a system with an age similar to that of Kepler-78 
to have a planet inside 0.01 au. but it is more likely for systems older than 
1 Gyr, than for those with ages below 1 Gyr. 


1 Gyr, Figure [6] suggests that maybe as many as 4 out of 10000 
could survive inside a = 0.01 au. Repeating the calculation above 
suggests that maybe 4 — 5 stars with ages similar to that of Kepler- 
78 could host such a close-in planet. Again, given the high transit 
probability for such a close-in system, detecting a planet such as 
Kepler-78b becomes possible. Our results therefore suggest that it 
is possible for this process to have produced a planet like Kepler- 
78b. Of course, if Kepler-78 is closer in age to 600 Myr, than to 
900 Myr, Figure[6]suggests that it would become less likely. 

4.3 Perturber properties 

The results above suggest that it is possible for an outer perturber to 
drive a Kepler-78b-like planet into a close-in orbit within the age of 
Kepler-78. To see how the properties of the outer body influences 
the inner planet, we show - in Figure [7]- how the final semimajor 
axis of the planet depends on the mass of the outer body. Again, we 
only shows results from the simulation with Q p = 5000 and Q’ s = 
5 x 10 6 . Figure[7]suggests that there isn’t a particularly strong mass 
dependence, consistent with our simulations e ssentially being in 
the test particle regime dLithwick & Naozl201ll) . However, Kepler- 
78, which has an apparent magnitude of m v = 12, is not known 
to host a stellar companion. Since Kepler is sensitive down to an 
apparent magnitude of m v = 14, that would suggest that if there 
is an undetected companion it would need to have a mass less that 
about 0.5 Mq. 

Figure [8] shows how the final semimajor axis of the planet, a, 
depends on the orbital properties of the outer body. The top panel 
shows that it is more likely that the planet will end up close to the 
parent star, if the outer body is in a relatively close orbit ( a 0 ~< 
100 au). Kepler’s has a 4” pixel size and so Figure[7]does suggest 
that a sufficiently faint, non-variable companion - that could have 
pertubed a planet into a Kepler-78b-like orbit - could indeed have 
avoided detection. The bottom panel of Figure [8] shows how the 
final semimajor axis of the planet depends on the on the eccentricity 
(e D ) of the outer body. It indicates that close-in orbits are a little 
more likely when the companion has a high eccentricity (e > 0.4) 
but are still possible for those with smaller eccentricities. 
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Figure 7. Figure showing the planet’s final semi-major axis plotted against 
the mass of the outer companion. There appears to be little dependence on 
companion mass. Kepler would probably have detected a companion with 
a mass in excess of ~ 0.5 Mq, but this figure does show that a low-mass 
companion (M < 0.5 Mq) could indeed have produced a system like 
Kepler-78b system. 


Figure 9. Figure showing the obliquity of the inner system. The inner sys¬ 
tem starts with the orbital angular momentum aligned with the spins of the 
central star and planet. The perturbation from the outer companion can, 
however, cause the inclination of the inner orbit to oscillate and those sys¬ 
tem that are tidally circularised can end up with a range of obliquities. 



0.01 0.10 1.00 
a (AU) 


4.4 The pile-up inside 0.1 au 

Figures [3] and \5 \ show a pile-up of planets inside a = 0.1 au, peak¬ 
ing at ~ 0.02 au. In our simulations, between 350 and 850 (be¬ 
tween 3.5% and 8.5% of the full sample of 10000) had final semi¬ 
major axes inside a = 0.1 au (and had not reached their Roche 
limit). If we assume that 40% of the Kepler sample could have a 
binary companion (either primordial or through an exchange inter¬ 
action) - and that most Sun-like star s have terrestrial-mass, plane¬ 
tary companions dCassan et al . 20121 : [Greaves & Ricell201 111 - then 
our results suggest that as much as 3 % of the Kepler sample might 
have planets that have been perturbed into close-in orbits, with a 
distribution that pe aks at about 0.02 au. This is intriguingly similar 
to the suggestion in lSanchis-Oieda et all d2014l) that about 1 in 200 
Kepler stars hosts a planet with a period of 1 day or less. 



0.01 0.10 1.00 
a (AU) 


Figure 8. Figure showing how the orbital properties of the outer companion 
influences the final semimajor axis of the inner planet. The top panel shows 
that outer companions with smaller semi-major axes (a„) are more likely to 
drive the planet to within a = 0.01 au. The bottom panel shows that outer 
companions with large eccentricities are more likely to perturb inner planets 
into very-close-in orbits, but that it is still possible for outer perturbers with 
low eccentricities. 


4.5 Obliquity 

An interesting aspect of the Kozai-Lidov process is that it can per¬ 
turb a planet into an orbit that is inclined with respect to its initial 
plane and, hence, inclined with respect to the spin of the hos t star 
i iwu. Murray & Ramsahai|[2007l : iFabrvckv & Tremainell2007l) . We 
now have a number of close-in, ‘hot' Jupiters that are o n orbits in- 
clined with respec t to the spin of the host star dHebrard et al.l2008l : 
ITriaud et al.]|20ld ~) and these are thought to be a consequence of 
Kozai-Lidov cycles. Figure [9] shows the final angle (obliquity) be¬ 
tween the angular momentum vector of the inner planet’s orbit and 
the spin of the central star, and shows that a wide range of obliqui¬ 
ties are possible. All the systems initially have obliquities of zero 
(the angular momentum of the inner orbit is aligned with the spins 
of the parent star and planet) and Figure[9]shows that those that are 
perturbed into an inner orbit can then be tidally circularised with 
a large range of obliquities, consistent with other sim ilar studies 
dFabiwckv & Tremaineil2007l : Naoz^ Far r & Rasiolbo 1 2|) . Using the 
Rossiter-McLaughlin method to determine such a mis-alignment 
(e.g. IQueloz et al.ll2000h is probably not possible for such a low- 
mass planet, but it may be possible to do so using spot-cr ossing 
(De_sert et al. 201 lj) or astro-seismology IChaplin et al J|20 l3i) . 
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5 DISCUSSION AND CONCLUSIONS 

We have considered, here, if systems like Kepler-78b (an Earth-like 
exoplanet with a very close-in orbit) could be due to a perturbation 
from an outer companion on an, initially, inclined orbit. To do this, 
we consider a system in which the star and planet have the sam e 
masses and radii as in the Kepler-78 syste m dSanchis-Oieda et all 
l2013l : IPepe et al.l2013l : lHoward et alj2013l) . but in which the planet 
initially has an almost circular orbit with a semi-major axis between 
0.5 and 2 au. The system is also assumed to have an outer compan¬ 
ion with a semimajor axis between 40 and 20000 au, with an orbital 
eccentricity that can be as high as e = 1 (but constrained by sta¬ 
bility criteria) and that may be inclined with respect to the plane of 
the inner orbit. 

We ran a suite of Monte Carlo simulations in which we ran¬ 
domly select the inner and outer systems semi-major axes, the ec¬ 
centricity of the outer system, the mass of the outer companion and 
the mutual inclination of the two orbits. We ran two sets of simula¬ 
tions, one where each system was evolved for t = 800 Myr, similar 
to the expected age of the Kepler-78 system, and the other where 
the age of the system was randomly selected to be between 500 
Myr and 2 Gyr. Our basic results are that: 

• it is possible for a planet to be perturbed into an orbit similar 
to that of Kepler-78b around a star with an age (600 — 900 Myr) 
similar to that of Kepler-78. Out of a sample of 10000, between 4 
and 10 survive inside a = 0.01 au. 

• if we consider a broader age range, the likely binarity of the 
Kepler sample, and the size of the Kepler sample, our results sug¬ 
gest that as many as 20 of the Kepler targets with ages less than 2 
Gyr could host a Kepler-78b-like planet. Additionally, we find that 
a system with an age similar to that of Kepler-78 could indeed have 
been found to host a Kepler-78b-like planet. 

• a planet such as Kepler-78b will, quite quickly, reach its Roche 
limit and be tidally destroyed. Our results suggest that such a planet 
would only survive inside a = 0.01 au for a few hundred Myrs, 
at most. In the simulations here, typically in excess of 100, but 
no more than 300, (out of 10000) reached their Roche limit and 
were assumed to be destroyed. This appears consistent with other 
work that has also suggested that this process could lead to the tidal 
destruction of perturbed planets l lNaoz. Farr & Rasioll2012l) . 

• given Kepler’s 4” pixel size and magnitude limit, it is possible 
that a faint, non-variable companion that could drive Kozai-Lidov 
cycles may have gone undetected. That the companion appears to 
need to be inside 100 au, means that it may be possible to detect 
the resulting radial velocity drift. 

• If a planet such as Kepler-78b were perturbed into its current 
orbit through Kozai-Lidov cycles, we might expect the star’s rota¬ 
tion axis to be misaligned with respect to the planet’s orbit. Mea¬ 
suring the star’s obliquit y is quite difficult, but could be possible 
using spot-crossing dDesert et al.ll201 lh or using astro-seismology 
dChaplin et alj2013l). 

• even though it appears possible that a system such as Kepler- 
78b could form in this way, it appears to be more likely for system 
older than 1 Gyr, than for systems younger than 1 Gyr. 

Our basic results, therefore, suggest that such a process could 
operate, but there are some caveats. Although it is possible for a 
system with an age similar to that of Kepler-78 to host a planet 
like Kepler-78b, the numbers are small (we may expect the Ke¬ 
pler sample to host only a few such planets). Additionally it seems 
that it would have been more likely to have found such a planet 
in a slightly older system. These results, therefore, suggest that 


Kozai-Lidov cycles could have played a role in the evolution of 
Kepler-78b, but don’t rule out that there co uld be a n alt ernative ex¬ 
planation, such as planet-planet scattering dRasio etahll 19961) .This, 
however, may suffer from the similar issues, since the dominant 
constraint - given the relatively low age of Kepler-78 - is the tidal 
evolution timescale. This constra int would proba bly also apply to 
the tidal downsizing hypothesis dNavakshinll20ld ). in which a mas¬ 
sive gas-giant planet formed in the outer parts of the system, mi¬ 
grates rapidly inwards and loses masses via tidal stri pping . That , 
of course, leaves the possibility that disc migration dWardl 119971) 
moved this planet into a very close orbit which has since evolved, 
through tidal interactions with the host star, into the orbit it inhabits 
today. Again, this would also involve tidal evolution once the disc 
has dispersed and so may also have a similar timescale issue, unless 
disc migration can place the planet sufficiently close to the parent 
star so that it can the n tidally evolve to where it is today. 

Recent work bv ISanchis-Oieda et aid d2014l) suggest that about 
1 out of every 200 stars hosts an ultra short period (USP) planet (pe¬ 
riod of 1 day or less). Although we’ve focussed on Kepler-78b here 
and found that few of our simulated systems have final periods as 
short as Kepler-78b (8.5 hours), many more have periods of 1 day 
or less. The exact number depends on the chosen parameters, but it 
varies from ~ 100 to just over 300 (from a total sample of 10000). 
Given that not all stars are binaries, th is is intriguingly similar to 
the result in lSanchis-Oieda et al.l | |2014|) . Similarly, our results sug¬ 
gest that such a process should lead to a pile -up of pla nets with a 
peak a t about 0.02 au, again consistent with ISanchis-Oieda et al.l 
d2014l) who find that the occurence rate rises with period from 0.2 
to 1 day. It may, however, be difficult to distinguish a pile-up due to 
Kozai-Lidov c ycles from what is expected from scattering in multi¬ 
planet svstems fSchlaufman, Lin & Ida! d20ld) . ISanchis-Oieda et all 
d2014l) do, however, suggest that almost all USPs have companion 
planets with period P < 50 days, which may provide a constraint 
on the formation process for these USPs. 

We should also acknowledge the possibility that our assump¬ 
tions do not properly represent the possible initial conditions such 
a system could have. The initial distribution of the planet in semi¬ 
major axis space may be different to what we’ve assumed and the 
orbital properties of the outer perturber may also be different. Sim¬ 
ilarly, the tidal properties of the parent star and planet may dif¬ 
fer from what we’ve assumed. However, we should at least ac¬ 
knowledge that even though our results suggest that Kozai-Lidov 
cycles will rarely produce a planet with properties similar to that of 
Kepler-78b, Kepler-78b is itself rare. In that sense our results could 
be seen as somewhat consistent with our knowledge of such plan¬ 
ets, but that - alone - doesn’t allow us to determine if it is likely that 
such a process did indeed play a role in the evolution of Kepler-78b. 
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